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Chemical reactions inside cells occur in compartment volumes in the range of atto- to femtolitres. 
Physiological concentrations realized in such small volumes imply low copy numbers of interacting 
molecules with the consequence of considerable fluctuations in the concentrations. In contrast, rate 
equation models are based on the implicit assumption of infinitely large numbers of interacting 
molecules, or equivalently, that reactions occur in infinite volumes at constant macroscopic concen- 
trations. In this article we compute the finite-volume corrections (or equivalently the finite copy 
number corrections) to the solutions of the rate equations for chemical reaction networks composed 
of arbitrarily large numbers of enzyme-catalyzed reactions which are confined inside a small sub- 
cellular compartment. This is achieved by applying a mesoscopic version of the quasi-steady state 
assumption to the exact Fokker-Planck equation associated with the Poisson Representation of the 
chemical master equation. The procedure yields impressively simple and compact expressions for 
the finite-volume corrections. We prove that the predictions of the rate equations will always under- 
estimate the actual steady-state substrate concentrations for an enzyme-reaction network confined 
in a small volume. In particular we show that the finite- volume corrections increase with decreasing 
sub-cellular volume, decreasing Michaelis-Menten constants and increasing enzyme saturation. The 
magnitude of the corrections depends sensitively on the topology of the network. The predictions 
of the theory are shown to be in excellent agreement with stochastic simulations for two types of 
networks typically associated with protein methylation and metabolism. 



I. INTRODUCTION 

Recent years have seen a distinctive surge in the formu- 
lation and application of stochastic models of biochemical 
reaction kinetics. This trend has resulted from a deeper, 
ongoing appreciation of the conditions characteristic of 
the intracellular environment^ and of their dissimilarity 
from in vitro conditions. Typical in vivo concentrations 
are in the range of nanomolar to millimolar; such con- 
centrations realized in a macroscopic volume imply very 
large copy numbers of interacting molecules whereas the 
same concentrations in the small volume of a cell fre- 
quently implies copy numbers ranging from few tens to 
at most few thousands (for a detailed experimental pro- 
tein abundance study see, for example, Ref. 2). Reaction 
kinetics is inherently a stochastic process;'^ this noisiness 
is not apparent in macroscopic conditions due to an im- 
plicit averaging over a very large number of molecules 
but cannot be overlooked when we are studying the ki- 
netics of a system in which the copy number of at least 
one species is small. This is frequently the case of intra- 
cellular kinetics. 

The introduction of the stochastic simulation algo- 
rithm by Gillespie^ has popularized the numerical study 
of stochastic reaction kinetics. However to-date the an- 
alytical study of the properties of such systems has re- 
ceived comparatively very little attention principally be- 
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cause the mathematical formalism of stochastic kinetics 
[i.e., chemical master equations, (CMEs)] is very differ- 
ent than that of deterministic kinetics [i.e., rate equations 
(REs) which are based on ordinary differential equations] 
and is less amenable to analysis. This problem is aug- 
mented by the fact that many biological networks of in- 
terest are considerably large. 

One of the main analytical methods for systematically 
exploring the stochastic properties of these networks has 
been the linear-noise approximation.^'^ The advantage of 
this method is the relative ease with which one can com- 
pute the magnitude of intrinsic noise (i.e., coefficients of 
variation and Fano factors). The major drawback is that 
the linear-noise approximation gives only meaningful re- 
sults provided the copy number of molecules is not small 
or to be more precise it is correct in the limit of infinitely 
large reaction volumes, i.e., the same limit in which the 
REs are valid. Since intracellular reactions occur in the 
opposite limit of small volumes, it is highly desirable to 
calculate the finite-volume corrections to the concentra- 
tions and moments of intrinsic noise. 

Developing a theory of finite-volume corrections 
presents a considerable analytical challenge. One way 
to obtain the latter is via the system-size expansion of 
the CME.^ The linear-noise approximation comes about 
by evaluating the first term (of order V'^ where V is the 
reaction volume) in this expansion. The next term is pro- 
portional to V^^i/2 and hence consideration of this term 
will necessarily give a finite-volume correction. Grima '® 
calculated the first such corrections for the mean concen- 
trations of species involved in single-substrate enzyme re- 
actions and recently also for a general chemical reaction 
network of arbitrary complexity.^ The general result is 



that the kinetics of a rcaetion pathway confined in a vol- 
ume V can be described by the usual REs plus new terms 
which are proportional to V~^. These equations are re- 
ferred to as effective mesoscopic rate equations (EMREs). 
The differences between the solution of EMREs and the 
corresponding REs for finite volumes stems from a cou- 
pling between the mean concentrations and the fiuctua- 
tions about them. EMREs can be explicitly solved for 
pathways characterized by a handful of chemical species, 
but otherwise one has to resort to numerical solution. In 
the latter cases, one obtains a solution at the expense of 
losing the insight which typically comes from analytical 
results. 

In this article we develop an alternative, powerful 
method of calculating finite-volume corrections to the 
solutions of the REs. Note that in the context of this 
article, finite- volume corrections exclusively refer to cor- 
rections to the mean concentrations not to the moments 
of intrinsic noise. The method is based on the Pois- 
son representation of the CME rather than the system- 
size expansion used in the derivation of EMREs. This 
method unlike the system-size expansion has been ap- 
plied to study systems of biochemical or biological rele- 
vance in only a handful of cases (see, for example, Refs. 10 
and 11) but as we shall show it is a tool with great 
potential for this field. We focus on chemical reaction 
networks which are composed of enzyme-catalyzed re- 
actions, a commonly encountered case in intracellular 
biochemistry. -"^^ We show that when the timescales of 
complex and substrate fluctuations are well-separated, 
it is possible to obtain explicit and impressively simple 
equations for the finite-volume corrections. For simple 
reactions these corrections are shown to be the same as 
given by the EMRE. The distinct advantage of the new 
method over the EMRE is that it provides analytically 
simple results even for complex networks with hundreds 
or thousands of species. This is inherently possible be- 
cause of the large reduction in the effective dimensional- 
ity of the CME when timescales are well separated. 

The paper is organized as follows. In Section II, we 
derive the Poisson Representation for a general enzyme- 
reaction network and use the resulting Fokker-Planck 
equation (FPE) to obtain an exact Liouville equation 
encoding all information about deviations from the de- 
terministic solution of the REs. In Section III we show 
that in the limit of well-separated timescales of complex 
and substrate species, the Liouville equation simplifies 
to a compact approximate form. This is used in Section 
IV to compute explicit expressions for the finite-volume 
corrections of a general network. In the latter section we 
show that the corrections for a simple Michaelis-Menten 
type reaction agree with those previously derived using 
the EMRE formalism. More importantly we apply the 
theoretical results to two common types of large-scale 
networks and confirm the predictions using simulations. 
Wc finish by a discussion in section V. 



II. THE POISSON REPRESENTATION FOR 
THE ENZYME REACTION NETWORK 

In this section we use the Poisson Representation to 
derive a general FPE for an enzyme reaction network. 
The latter while being exactly equivalent to the CME is 
much more amenable to analysis and hence is a very con- 
venient starting point for detailed calculation purposes. 

We consider a generic type of enzyme network com- 
posed of two major types of chemical processes: (i) the 
input of a substrate species Aq into a subcellular com- 
partment; (ii) the transformation of Aq into some final 
product Atv via N consecutive enzyme-catalyzed reac- 
tions of the type 



feL, 



(1) 



where Ai, Ci and Ei denote the i^^ substrate, complex 
and enzyme species, respectively; the index i takes values 
from to A^ — 1 and the fc's denote the relevant macro- 
scopic rate constants. Note that we have assumed here 
that the bimolecular reaction rate, fci , is the same for all 
substrates. Note also that in such types of networks there 
are N distinct substrate species, an equal number of dis- 
tinct complex species and a number of enzyme species 
varying between 1 and TV. The freedom in choosing the 
number of enzyme species comes from the fact that an 
enzyme can generally bind to more than one type of sub- 
strate. 

If we assume that we have well-mixed conditions inside 
the compartment then the instantaneous description of 
the state of a chemical system at time t is simply given 
by the vector of the absolute number of molecules of each 
species, n = {{nAi}, {^Ei}, {nci})- Since the mesoscopic 
kinetics are stochastic, a full description of the system is 
necessarily probabilistic and is achieved by defining the 
probability density function, P — P{n, t) and its time- 
evolution equation, which is commonly referred to as the 
CME (Refs. 5 and 13) 



N-l 

dtP = ^ n,p + EP, 

i=0 

k^ 
Ili = -rr(0A.0Q - l)nA,nEi 

+ Fi(ec.e^i-l)nc. 
+ ^2(0c.e^i^^-l)nc,, 
S = h^V{Q^^ - 1), 



(2) 



(3) 
(4) 



where V is the compartment volume, S is the contribu- 
tion due to the input of substrate species Ao into the 
system at a rate km while Hi describes the «*'' catalytic 
reaction step, Eq. (1), in the reaction network. The 
CME is compactly expressed using van Kampen's step 
operators defined as Q'^^g^nxi) = g{nxi ± !)■ 

Substantially, the CME depends only on the set of vari- 
ables {riAil and {nci} since an enzyme molecule can be 



either in the free state or in the complexed state and 
hence the variables {n^;.} are redundant. We can express 
this conservation law by writing ue^ — E^ — J^j GijUCj 
where the matrix Gy is defined by construction to be 



G, 



(ei)j 



f , if enzyme i binds substrate j 
0, otherwise. 



(5) 



The A'^-dimensional vector, Cj, is associated with the en- 
zyme binding substrate in the i-th catalytic reaction. Its 
j-th entry is chosen to be equal to one if the enzyme can 
form a complex with the j'-th substrate and zero other- 
wise. Hence the connectivity of the network is explicitly 
encoded in the form of these vectors. 

We define the moment generating function, param- 
eterized by the vector of continuous variables, z = 
{{^aAAzc,}), as 



Gw-En 






(6) 



Muhiplying the CME, Eq. (2), by Hi ^Z' ^Ci ' perform- 
ing the summation over all values of the variables n and 
expressing the resulting equation in terms of the moment 
generating function G{z)^ we obtain the moment gener- 
ating function equation 



5« = K,V{zAo - 1), 



TZf ^ {zc. - ZA^) (£^d^^ - Kljdt' 



+ Kl{zA^^,-ZA^)d^^, 



(7) 
(8) 

(9) 



where £i abbreviates [E^] — V^^Tl^Cid};; and [Elj] = 
E^/V is the total enzyme concentration associated with 
the enzyme binding substrate in the i*'' catalytic reaction 
step. Furthermore we set Kl ~ k^_-^/ki, K2 = k\/ki, 
K\i = K\ + K2 (commonly referred to as the Michaclis- 
Menten constant) and 9^ = d/dzx- Note that in Eq. (9) 
all derivatives are to the right. 

We proceed by making use of Gardiner's Poisson Rep- 
resentation. At the heart of this method is the assump- 
tion that the probability density function P(n,i) can be 
expanded as a superposition of multivariate uncorrclated 
Poissons^'^'^^ 



P{n,t) = / da Jl 



,-aA,V 



{o^A^vy 



"Ai! 






fia,t), (10) 



where the function f{ck,t) is usually referred to as the 
quasi-probability density function. The vector a is de- 
fined to be ({a^i;}, {aci})- It has been shown that this 
superposition always exists if the range of a is extended 
to the complex plane by analytic continuation of the Pois- 
son kernel. We have explicitly introduced V in our su- 
perposition definition, which is not customarily done in 



the Poisson Representation. In doing so, we obtain the 
representation in intensive variables, i.e. in units of con- 
centrations, as those encountered in the theory of rate 
equations. The above expansion is equivalent to writing 
the moment generating function as 



G{z,t) = I daf{a,t)e 



i:{zAi-l)aA,V ^j:(zc^~l)aCiV 



(11) 



It follows from Eqs. (7) and (11) (see Appendix A for 
details) 



-dtf = in + 5)/, 

o = kind A, 

7^ = fcl E^^' 

4=0 



(12) 
(13) 



^ Kl{d^' - d\)ac., (14) 

where we have utilized the notation v = V^^^^, dx = 
d/daxi- Note that aA„ = cep is the variable for the 
product formed after N catalytic steps. Note also that 
£i is not a constant, but is given by the operator 



£^ = [Ei.]-J2G.,il-,.^^c,)ac,, 



(15) 



which is essentially the Poisson representation of the con- 
servation of total enzyme molecules. The above equa- 
tion generally differs from the corresponding determinis- 
tic conservation law in terms of concentrations; the latter 
is described only by its average, while the former exhibits 
finite volume corrections due to the finite copy number 
of enzyme molecules. Given Eq. (15), we see that the 
Poisson representation, Eq. (12), yields a Fokker-Planck 
equation in terms of substrate and complex variables. 
Note also that for the case of A^ = 1 we obtain the repre- 
sentation for the single-substrate single-enzyme reaction, 
which is usually referred to as the Michaelis-Menten re- 
action. 

This completes the derivation of the Poisson represen- 
tation of our general enzyme-reaction network. Note that 
this is not the same FPE as that which arises from the 
system-size expansion method of van Kampen.^ In the 
latter case, the FPE is an approximation to the CME 
in the limit of large volumes whereas the FPE obtained 
from the Poisson Representation is exactly equivalent to 
the CME. 

A further boon of the Poisson representation is that 
once we have calculated the moments of the continuous 
variables, axi, using the FPE, we can very easily find the 
corresponding moments of the copy number of molecules, 
nxi, using the following simple relationships:^^ 



V 

{nx,nxj) 

y2 



("X,), 



(16) 



{ax,ax,) + —5ij{ax,), (17) 



where the angled brackets imply the statistical average: 
on the left hand side of the above equations these are 
given by (..) = J dn .. P{n,t) while those on the right 
hand side imply (..) = J da .. f{a,t). 



A. The mesoscopic equation 

The kinetics becomes deterministic in the macroscopic 
limit of infinitely large volumes. This can be easily ver- 
ified by noting that all second-order derivatives in the 
FPE are multiplied by a factor proportional to the inverse 
of the volume. To compute the finite- volume corrections 
we will need to separate the mesoscopic and macroscopic 
evolution equations. We now show that this can be done 
by applying a suitable change of variables to the Pois- 
son Representation. The macroscopic corresponds to a 
shot noise contribution, which agrees on average with the 
mean-field result i.e. with the solution of the REs for the 
substrate-enzyme network. The mesoscopic contribution 
refiects the non-equilibrium properties of the network due 
to the bimolecular character of the substrate-enzyme in- 
teraction, and depends parametrically on the mean-field 
expectation. 

We express the deviation from the deterministic path 
by the following change of variables: 

{aA^,ac,) -> ([^^](i) + i^eA., [C,]{t) + lyec,) (18) 

We shall refer to ex; as the mesoscopic correction to 
the deterministic, macroscopic concentration [Xi](t) of 
species Xi. Note that the above equation has the same 
apparent form as the van Kampen (VK) ansatz,^'^ at 
the heart of the system-size expansion, but the context 
of the application is completely different. The VK ansatz 
is applied to the integer number of particles in the CME 
leading to an infinite series in powers of the inverse square 
root of the volume, which has to be truncated; the first 
term of this expansion (the one proportional to V'^ ) gives 
a linear FPE which is an approximation to the CME. In 
our case the change of variables, Eq. (18), is applied on 
the FPE arising from the Poisson representation which 
leads to a finite series and allows for exact analytical 
treatment; as we shall show now, this divides the exact 
FPE into a macroscopic term and a term which captures 
all deviations from the macroscopic. 

The transformation applied to the FPE transforms the 
time derivative into 
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.f{[X]+ue,t) 

V,f{[X] + ue,t), (19) 



where a/at\x denotes taking the derivative with x held 
constant. It follows from Eq. (18) that de/dijo, = 
—i'~^d[X]/dt. Finally by expressing the right hand side 
of Eq. (12) in terms of the new variables e and equating 



the result to Eq. (19), we find that the FPE takes the 
form 

-—g{e) ^kiCg[e) + t^ M ki'Rnia.ao + 5macro 

^M v„-ffl^vj,(.). ,20, 



at 



at 



where the relevant operators are 



5, 



macro — "'in 



— kw 



a 



de, 



(21) 
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dec. 



i 

+ Y,{Kl[C,]~m[A,])^, (22) 

i ^ 

^ = 5Z,{(^c - d\){[E,]eA, + [A,]5£, + ueA,5£^ 

- Al,6c.) + KHa^^' - a\)ec,). (23) 

Note that the new probability density function necessar- 
ily satisfies ,9(e)de = f{a)da.. Note also that the nota- 
tion 9 Y now denotes the derivative a/aexi ■ The quantity 
[Ei] ~ [E'!p] — J2n Gij[Cj] is the macroscopic concentra- 
tion of the free enzyme species associated with binding 
substrate in the i*'* catalytic reaction step. The operator 
6£. = -m] - £,)/u = - E, G,j{ec, - 9^[C,] - i^^^ec,) 
is the contribution of the enzyme-operator Eq. (15). 

Note that the resulting form of the FPE is clearly di- 
vided into two parts. In the macroscopic limit the terms 
proportional to u~^ dominate and their sum must equate 
to zero - this leads to the macroscopic equations. It also 
then follows that the mesoscopic equation is simply given 
by 



arg{e,t) = -Cg{e,t), 



(24) 



where we have renormalized time to r = kit. Note also 
that due to the Poissonian nature of the substrate in- 
put process, it only contributes to the macroscopic part, 
a feature which is unique to the Poisson representation. 
We emphasize that up till this point, we have made no 
approximations and hence the resulting mesoscopic equa- 
tion is exact. 



III. ADIABATIC ELIMINATION OF THE 
COMPLEX SPECIES VARIABLES 

In this section we show how to rigorously eliminate the 
fast variables from our description. This will be done in 
two steps: on the macroscopic contribution of the FPE 
and on the mesoscopic contribution including terms up 
to order i/ (i.e., finite volume corrections). As we shall see 
later on, the reduced mesoscopic description does depend 
on the reduced macroscopic description and hence the 
need to treat the latter first. 



A. The reduced macroscopic equations 

The macroscopic equations arc obtained from Eq. (20) 
by taking the limit i^ — > 0, leading to 



dt 

d[Q] 

dt 

d[AN] 

dt 



k,[E,][A,]-{ki + ki,)[a], 

^2 " [Cn-i], 



[Ei] = [E'^]-22 ^iJ l^jl- 



(25) 



These agree exactly with those that can be obtained from 
the RE approach. The elimination of the complex species 
from the macroscopic equations is a well-known proce- 
dure commonly referred to as the quasi-steady state ap- 
proximation (QSSA).^^'^^ Briefly speaking the approx- 
imation is tantamount to assuming that the complex 
equilibrates on a much shorter timescale than the sub- 
strate. This is implemented by imposing the approxima- 
tion d[Ci]/dt = on the macroscopic equations. The re- 
sulting reduced equations (commonly referred to as the 
Briggs-Haldane equations) are then given by replacing 
[C'i] in the full time-evolution equations for the substrate 
concentrations by 



[C.l- 



m[A, 



K 



(26) 



M 



B. The reduced mesoscopic equation 

Now we are interested in deriving the reduced meso- 
scopic equation corresponding to the reduced macro- 
scopic equations that we just considered. Time-scale sep- 
aration on the mesoscopic scale is non-trivial because of 
the inherent correlations between the mesoscopic fluctu- 
ations of the various species. Our presentation shall be 
as follows. First we shall show that the mesoscopic Liou- 
villian, Eq. (23), can be generally cast into an asymptotic 
form of the interaction representation which is typically 
encountered in the theory of adiabatic elimination of fast 
fluctuating variables. ^^^^^ The latter yields a particularly 
simple result for the reduced mesoscopic equation. 

We will now show that the Liouvillian can be rewritten 
in the general form 



£(7) = 7/:i + 7'/'/:2 + /:3 



(27) 



where 7"^ is the characteristic fast timescale of the com- 
plex fluctuations which will be specified later on. We 
now proceed to derive the operators £1, £2, £3 for the 
general enzyme-reaction network under study. We start 
by grouping all terms containing only the pair of complex 
variables ({eciji {f^cD i^^*-" '^i' terms concerning solely 



the substrate ones ({cAi}, {^^i}) into £3, while treating 
the remaining terms as interaction £2. Thus we have 



A=^£«, C, = Y.ct\ C.^Y.^ 



(28) 



<i) 



C'l> =d^c{[M5e^^K\,ecd, 



<i) 



<i) 



-d\{[A,]6£,-K\,ec,) + K'2{d_ 



i+l 



(29) 

(30) 

(31) 



It is instructive to rescale all complex variables by their 
characteristic timescale 7 



z^ = l"^^c. 



^Ai- 



The operator £3 is trivially obtained 



£. 



('■) 



&l\Ei 



(32) 



(33) 



Note that dl. denotes the derivative d/dxi whereas d\ 
stands for the derivative d/dexi- Next we observe that 
the enzyme operator transforms as 



7 



''H£, ^-J2 Gvi^, idiiC,] v^^l-'dlz,). (34) 



,(j) 



(0 



Plugging this into Eq. (29) and putting C\ — > 7£^ 
we find that the last term in (34) can be asymptotically 
neglected. Hence the dominant contribution in the limit 
of large 7 is given by 



where we have utilized the abbreviations 



(35) 



(36) 



Thus the asymptotic form of the complex fluctuations 
is described by an Ornstein-Uhlenbeck process, centered 
on the deterministic expectation Eq. (26). It is clear, 
by virtue of the fluctuation-dissipation theorem, that the 
matrix J (constant matrices, i.e., those independent of 
the e^. and ed variables and of the associated partial 
derivatives, are underlined throughout the rest of the ar- 
ticle) must correspond to the Jacobian of the complex 
species while the matrix D determines the strength of 
the fluctuations in the Poisson representation. Note that 
throughout the article by Jacobian of a species we mean 
the negative of the Jacobian matrix as obtained from the 
macroscopic rate equations of that species with the time 
scaling r = kit. 

The characteristic timescale 7"^ can be inferred from 
the explicit form of the Jacobian. We choose 7 = tr(J) 
giving the relaxation rate of the complex vector z. Con- 
sequently M is flnite, as required by the stability of £1, 
and is given by the Jacobian of unit trace. Equally we 



observe that the Jacobian of the substrate vector x is 
given by the diagonal matrix E — diag([£'i]), such that 
its relaxation rate is proportional to its trace. Thus under 
the condition where 7 3> tr(E), we can assume that the 
timescales of substrate and complexes are well-separated. 
The result can be interpreted as follows: Ci represents 
the fast noncquilibrium fluctuations of the complex and 
£3 reflects the slow Poissonian decay of the substrate 
fluctuations. 

Finally we turn our attention to the interaction, we put 
^2 -^ 1^''^ 1^1 and rearrange to emphasize its dependence 
on the complex variables 

^ — 'j '^ — 'j 

^,, = [k, {E^\x, - 9^ A, - vS^^x^G,, [C,]) . (37) 

For notational convenience, we have considered only 
terms proportional to z and Vz, namely those which give 
non- vanishing contributions in the following. We empha- 
size that Eqs. (33), (37) together with (35) yield the con- 
tribution to the mesoscopic equation under the condition 
of well-separated substrate- and complex timescales. 

It is clear that the reduced substrate density function 
h{x) = J dzg{x, z), cannot simply be deduced from £3, 
but must retain a contribution from the inherent coupling 
to the rapid fluctuations of the complex species. How- 
ever, it is known^^ that in the limit of well-separated 
timescales, a particularly simple result can be obtained 
by projecting the evolution equation on the steady-state 
of the fast variables. The technique is the adiabatic elim- 
ination of the fast variables and can be accomplished by 
the use of the projector Pg = Tr{z) J dzg{x,z), with 
TT being the steady-state distribution of the fast vari- 
able vector z, determined by the steady-state condition 
£i7r = 0. 

The applicability of the method follows directly from 
the asymptotic form Eq. (27), which has been derived 
here explicitly for the large scale enzyme-reaction net- 
works under consideration, together with the conditions 
PjCs = C^P.PCi = CiP = 0,PC2P = 0,13 which can 
be easily verified from Eqs. (33), (35), and (37). The re- 
sult in the limit of large 7 has been formally derived by 
Gardincr^'^ and reads 

-dt{Pg) = lim PC{^)g = (£3 - PC2C^'C2){Pg). 

(38) 

Integration over z yields the reduced evolution equation 

~dth{x,t)^C'h{x,t), (39) 

C = lim £(7) = £3 - (£2£r'^2).. (40) 

7— foo 

The mesoscopic description is now reduced to that of the 
substrate only. The angled brackets (•)7r denote the trace 
/ dz • tt{z) taken over the rapid steady-state fluctuations 



of the complex species. We now substitute Eq. (37) into 
Eq. (40) 

C'^C,-Y.{x{zC^'z^).xn,, 



.^(x(^£rvr). 



^'>.y (41) 



The explicit form of the correlators involved has been 
derived in Appendix B using the steady-state condition. 
The result is 



{zC-^z'^)^ = S = M^^D^M"'^, 

-{zC^'S7^)^^M-\ 



(42) 



Hence using Eq. (41) together with Eq. (42), we get 
the final form for the reduced Liouvillian in a convenient 
matrix form, which yields the mesoscopic implementation 
of the QSSA 

£' = £3 + J2(x^X^)v + Y.ixM-'^'^h- (43) 



IV. FINITE VOLUME CORRECTIONS 

In this section we will use the reduced mesoscopic de- 
scription to compute the finite volume corrections to the 
macroscopic steady-state concentrations of all substrate 
species in the network. By definition, the concentration 
of species Aa according to the stochastic model is given 
by Pa = (na/V). Using Eq. (16) together with Eq. (18) 
it is straightforward to show that 



Pa = [Aa] +iy{Xa). 



(44) 



Hence the size of the finite-volume correction to the 
macroscopic concentration of species Aa is simply given 
by (xa)- The time-evolution equation for the latter can 
be computed directly from our reduced mesoscopic equa- 
tion, Eq. (39) together with Eq. (43), yielding 



dr{Xa) ~ / dxXaCh{x). 



(45) 



We evaluate the corrections associated with the sub- 
strate, a < iV, to leading order in x and v. Defining the 
transfer matrix T^ ^ = K2{Si^i^a — ^ia), one can show 
the following results: 



dxxC3h:^E{x), (46) 

dec a; ^(xM" V^)ij/i 

ij 

= (1 + TM"^)E(a;) + -MS', (47) 

7 

dx X J2ixSx^hh = - (M + T) 5, (48) 



where S is the vector obtained by summing over all rows 
of the matrix S, i.e., {S)i = ^ ■ Sf,. Summing up these 
equations, we obtain 



(49) 



dr{x) =T [M-''E{x)--S 



7 



which has the steady state solution (a;) = (iy/^)E~^MS. 
Note that the result is virtually independent of 7, if we 
rewrite 



(xa) = z/^(E-iDJ-^), 



(50) 



Hence it is clear that the finite volume corrections for 
substrate species are non-zero; this indeed signals the 
breakdown of the law of mass action on mesoscopic length 
scales or equivalently for low copy number of molecules. 
It can also be shown that the elements of the inverse 
Jacobian of the complex are given by 



^,^ 



1 



aiCri 



(51) 



where a = Aia.g{[Ai\ / K\j) and a^ is the i diagonal 
element of the latter matrix with value [Ai\/K\f (the 
reduced macroscopic substrate concentration of species 
i). Given the definition, Eq. (5), it is straightforward 
to verify from the above equation that the elements of 
the inverse Jacobian are always positive valued. Since 
both matrices D and E^^ are also positive it then fol- 
lows that the finite volume corrections to the substrate 
concentration, Eq. (50), are always positive. In other 
words, the predictions of the REs will always underes- 
timate the steady-state substrate concentrations for a 
substrate-enzyme network confined in a small volume. 

The equations for the finite volume corrections can 
be conveniently expressed in a form which is useful for 
obtaining insight about the physical origin of the cor- 
rections and also for their numerical computation. We 
write the matrix D ~ AGC, where A = diag([Ai]) and 
C = diag([Ci]). The Jacobian of the free enzyme is sim- 
ply the transpose of the Jacobian of the complex such 
that (J£j^)i = ^jJji^', this can be directly computed 
using Eq. (51). These two equations together with the 
steady-state condition, Eq. (26), written here as C = E a, 
the identity (E^^ G E)qj = {ea)j and Eqs. (44) and (50), 
allow us to write a simple, final equation for the steady- 
state substrate concentration as predicted by the stochas- 
tic model 



[Ao.] 



— 1-1- Ra, 



where 



1 



rCa — ^ ^^a—^ E ■ 



(52) 



(53) 



The parameter, R^, is proportional to the product of 
the relaxation time of the enzyme species and the corre- 
sponding reduced macroscopic substrate concentration, 



summed over all substrates which can bind to the same 
enzyme. It is also straightforward to derive expressions 
for the relative error made by the RE model, for the ab- 
solute differences between the reduced mesoscopic and 
macroscopic concentrations and for the absolute differ- 
ences between the sum of the reduced mesoscopic and 
macroscopic concentrations, respectively 



R. 



Pa - [Aa] Ro 



Pa 1 + Ra' 

Pa - [Ag] _ 



M/ 



E Pa - [^g] _ V^ n 



K" 
^M 



(54) 
(55) 
(56) 



In the following subsections, we will apply the general 
results developed so far, to three different cases which 
are commonly encountered in biology. We will verify the 
main theoretical predictions, i.e., Eqs. (54)-(56), by com- 
parison with detailed stochastic simulations. 



A. Single-substrate reaction with single-enzyme 
species 

The reaction where a single enzyme-species catalyzes 
only a single-substrate is the classical Michaelis-Menten 
reaction textbook example 



^0 , Eq + Aq ^f=^ Co 



Eo + Ai. 



(57) 



The enzyme only binds to a single substrate species and 
hence Goo ~ 1. It then follows that Eq. (53) evaluates to 



i?o — 



1 



flO 



K°jV{ao + l)' 



(58) 



where oq — [Ao]/K^j. Substituting the above in Eq. (52) 
gives an expression for the ratio of the mesoscopic and 
macroscopic substrate concentrations; this agrees exactly 
with that obtained by taking the limit of well-separated 
time scales (i.e., 7 = ([^0] + Km) > [Eo]) of the finite- 
volume correction previously calculated by Grima^ us- 
ing the system-size expansion including terms of order 
]/-i/2^ Note that the finite volume corrections are sig- 
nificant for large steady-state substrate concentrations, 
i.e., when the enzyme is working in saturated or near- 
saturated conditions. The predicted dependence of the 
relative error [i.e., Eq. (54) together with Eq. (58)] on the 
size of the finite volume was tested using stochastic sim- 
ulations (Fig. 1). Numerics and theory agree very well 
over at least 3 orders of magnitude of the compartment 
volume, or equivalently over the whole typical physiolog- 
ical range of enzyme copy numbers (from 1 to 1000). 
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FIG. 1: Single-substrate, single-enzyme reaction. Finite vol- 
ume scaling of the relative error in the substrate concentration 
as predicted by the RE model, i.e. -Rerror = (po — [^o])/po. 
The error stems from a finite-volume correction to the mean 
concentration due to intrinsic noise. The solid line shows 
the theoretical estimate as given by Eq. (54) together with 
Eq. (58) in the text. The data points are obtained from 
simulations carried out using Gillespie's stochastic simula- 
tion algorithm under steady-state conditions. The param- 
eters are: K'lf = 1/5 (fci = 5, fc^ = 0.5, A:°i = 0.5) 
and [E^] = 0.5. Timescale separation is guaranteed since 
7 — {[Ao] -\- KIj) = 10[i?ol and near-saturation conditions en- 
sue since fciii/(fc2 [£-§■]) = 0.8 (a value of 1 indicates complete 
saturation). The upper x-axis indicates the corresponding 
discrete number of enzyme molecules [-Ey]!^ used in the sim- 
ulation. Note that the agreement between theory and simula- 
tion is excellent even when there is just one enzyme molecule 
in the compartment. 



B. Multi-substrate network with single-enzyme 
species 

We now consider a network where different substrates 
compete for catalysis by a single-enzyme species 

-^ Aq, Eq + Aq n Co — ^ Eq + Ai, 
Eo + Ai^Ci^Eo + A2, 



Eq + An^i s C'n-I 



Eo + A 



N- 



(59) 



Such type of reactions are commonly associated with 
the methylations of proteins, DNA and RNA by the 
methyltransferase class of enzymes. A strictly unidirec- 
tional methylation occurs for instance in the nitrogen 
trimethylation of Lys^^. A single amino acid, embed- 
ded in a protein, can typically be methylated only a few 
times (e.g., three times in the case of Lys). But the 
mere large number of amino acids composing proteins 



means that the number of methylations can generally 
be very large. The modification of the protein from the 
unmethylated (or even partially methylated) form Aq to 
the fully methylated form An is processed by the enzyme 
within N catalytic steps. Partially methylated interme- 
diate states are then described by Ai,. . . ,An-i- This 
type of network has also been proposed as a candidate 
for explaining the kinetics of protein digestion. ^° 

The enzyme can now bind to all the substrate species 
and hence G^ = 1. It also follows that since we have one 
enzyme species then [Eip] = [E!^] = [Ex]- It can then be 
shown that Eq. (53) evaluates to 



i?, 



= -V — 



V^Ki,a + l 



(60) 



where a = Y^^ai is the sum of reduced macroscopic con- 
centrations. To obtain further insight into the above 
equation and to evaluate Eqs. (54)-(56) it is necessary to 
obtain expressions for the reduced macroscopic substrate 
concentrations. Using the REs, Eq. (25), and applying 
the QSSA to the complex concentration variables we find 
that the reduced macroscopic equations are 



d[AQ 
dt 



- k- - 7)° 



d[A, 
dt 



v'-^ - v' 



(61) 



where w' = fc2[^T] ai/(l + a), i.e., the catalytic velocity 
of step i. Imposing the steady-state condition on the 
macroscopic substrate concentration equations, Eq. (61), 
we find that the catalytic velocities of all steps must be 
the same and equal to the input rate, w* = fcjn. It follows 
that the reduced macroscopic substrate concentrations 
are a^ = r]i{l — r])~^ and their sum is a = rj{l — r])~^ 
where r]i is the dimensionless parameter defined as 



"^in 



Vi 



[E 



■T\kl 



(62) 



and ry = ^^ r]i . Note that rji is the ratio of the input rate 
and of the maximum rate at which the enzyme can cat- 
alyze substrate Ai into substrate Ai^i. Similarly rj is the 
ratio of the input rate and of the overall maximum rate 
at which the enzyme can catalyze the initial substrate Aq 
into the product Af^. Hence both rji and rj have a value 
between and 1. It follows also that 77 is a measure of 
enzyme saturation. 

Given the above results one can now evaluate Eqs. (55), 
(56), and (60) which leads to 



1 ^ ?7^- 

j ^M 



K 



M 



1 — ?/' 



y^ Pt - [^»] ^^ V 



K 



M 



1-r; 



(63) 

(64) 
(65) 



From Eqs. (52), (62), and Eq. (63) one can deduce that: 
(i) the concentration of each different substrate species is 
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FIG. 2: Scaling of the finite-volume corrections with network 
size for a multi- substrate network with single- enzyme species. 
Plot of the absolute difference between the sum of the reduced 
mesoscopic and macroscopic substrate concentrations versus 
the network size, A'^, for a homogeneous set of rate constants, 
fell = fc!^i, ^2 = fc§. Theory predicts that if the input rate is 
rescaled by the network size, fcin = kf^/N, then the absolute 
error in steady-state conditions [as given by Eq. (65)] should 
be independent of network size. Data points are obtained 
from stochastic simulations using Gillespie's algorithm (with 
a rescaling of kin as mentioned above) and solid lines are the 
theoretical estimates using Eq. (65). Note that in the y-axis 
label we have used p to mean J2i P» S'^d [At] to mean ^HJ^i]- 
Parameter values are found in the text. 



amplified (from the deterministic concentration) by the 
same factor; (ii) the deviation from the predictions of 
the deterministic model increases with decreasing vol- 
ume, decreasing Michaelis-Menten constants, increasing 
enzyme saturation and increasing network size N. The 
latter stems from the fact that increasing the number of 
substrates will increase competition for the single enzyme 
species meaning that the copy number of free enzyme 
molecules at any one time becomes smaller and hence 
leads to increased noise-induced effects. 

We have carried out extensive stochastic simulations 
using the Gillespie algorithm^ to test the accuracy of our 
predictions. The four tests are as follows: 

Test 1: Dependence of the finite-volume corrections 
with network size. We first impose homogeneous rate 
constants, i.e., a fixed set of rate constants for all en- 
zyme reactions in the network. It then follows that if we 
make the scaling ki^ — k?^/N, the quantity rj = ^^ rji will 
be independent of the network size. Consequently under 
this scaling, Ri , is predicted to be independent of the size 
and so is the sum of the absolute reduced concentrations, 
Eq. (65). Stochastic simulations show perfect agreement 
with this prediction (Fig. 2) which implicitly proves that 
the theory correctly predicts that finite volume correc- 
tions increase with network size. The relevant simulation 
parameters are: V = 500, ki = 1000, k'_^ = 0.5, fc| = 20, 
(KmV)-^ = 0.1, [Et]V = 5, fcP„ = 0.1, rj = 0.5 and 
7 == 10 tr(E). 

Test 2: Dependence of the relative error with en- 
zyme saturation. We fix the Michaelis-Menten constant. 



from Eqs. (54) and (63) that the relative error is simply 
given by R\„or — v/iv + KmV) i.e. the relative error 
increases with increasing enzyme saturation. The predic- 
tion is verified by stochastic simulation [Fig. 3(a)], where 
the rate constants were chosen so that k^ = /c^ for i even 
and ^2 = f^^k2 for i odd and the specific parameter val- 



V = 500, fci = 1000, N = 6, f 

== 10, 7 > 10 tr(E), {KmV)-' - 



= 0.5, 
0.05, 



ues used were: 

Kin ^ 1/0, 1^2 

[Et]V = 5. 

Test 3: Dependence of the relative error with species 
type. As previously mentioned, an inspection of Eq. (63), 
shows that the quantity Ri and hence the relative error 
Rlrror ^^^ predicted to be the same for all species in 
the network. This agrees with the results of stochastic 
simulation [inset of Fig. 3(b)] of a six substrate network 

where isTV = (1 + ^)^°/- ^2 = (1 + 0^2, ^-i = (l + O^-i, 
and the specific parameter values were: V = 1000, kin = 
1/60, k^ = 9, fc'^i = 1, fci = 1000, [Et]V = 5. 

Test 4-' Dependence of the finite-volume correction on 
the position in the network. Here we test the ability of 
the theory to predict the absolute difference between the 
reduced mesoscopic and macroscopic substrate concen- 
trations as predicted by the stochastic and deterministic 
models, i.e. Eq. (64), for each individual species in the 
network. For rate constants chosen as in Test 3, the dif- 
ferences are predicted to be proportional to (1 -f- i)~^, 
which is confirmed by simulations [Fig. 3(b)]. All pa- 
rameter values are exactly as in the previous test. 



C. Multi-substrate network with multi-enzyme 
species 

Finally we consider a sequential reaction network 
which is typically associated with metabolism^^ 



H Ao, Eo + Ao^Co^Eo + Ai, 
Ei+Ai^Ci -%E,+A2, 



En-1 + ^AT-l N Cn-1 - 



En-1+A, 



N- 



(66) 



This network is characterized by a different enzyme 
species for each separate catalytic step and thus Gij = 
dij. The finite- volume corrections can be easily derived 
from Eq. (53), since the Jacobian of the complex species 
is instantly diagonal and hence the result is 



Ri 



1 



K\,V (a. + 1) 



(67) 



K 



M 



Km, and the size of the network. It then follows 



It is useful to write the above expression in terms of the 
constants defining the network by using the REs (with 
the QSSA applied to the complex concentration vari- 
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Hence it follows that Eq. (67) reduces to the simple form 
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FIG. 3: Dependence of the deviations from the RE predic- 
tions with enzyme saturation and position in the network for 
a multi-substrate network with single-enzyme species. In (a) 
we plot the relative error as a function of rj which is a measure 
of enzyme saturation. In (b) we plot the absolute differences 
in the reduced mesoscopic and macroscopic substrate con- 
centrations as a function of the substrate species index and 
hence as a function of position in the six substrate network. 
The inset shows the relative error for all species. The data 
points are from stochastic simulation and the solid lines [in 
(a)] and the open circles [in (b)] are the theoretical predic- 
tions. The rate constants are heterogeneous in both cases 
but the Michaelis-Menten constants for each reaction in the 
network are the same for (a) while they vary according to po- 
sition in the network in (b). See the text for parameter values 
and for a detailed discussion. 



ables) 



where 



dt 



"-in i^M I 



d[A, 
dt 



^M 



^M' 



(68) 



'M 



k\[E'^\ ai/{l + Qi) is the catalytic velocity 
in this case. The connectivity of the network requires 
that all catalytic velocities must be the same and equal 
to the input rate, v\j = fcjn, from which it follows that 
the reduced macroscopic concentrations are 



ai 



m 



1 -?7i 



(69) 



R^ 



KmV 



m- 



(70) 



where r]i = ki^/k2[E^], a measure of saturation as before. 



It can be deduced from the above equation that the de- 
viation from the predictions of the deterministic model 
increases with decreasing volume, decreasing Michaelis- 
Menten constants and increasing enzyme saturation as 
for the previous network. However unlike the previ- 
ous case, for the metabolic network we find no depen- 
dence on the network size N and the factor Ri is species- 
specific implying that the concentration of each differ- 
ent substrate species is amplified (from the deterministic 
concentration) by a different factor. Hence the finite- 
volume corrections in this network are locally determined 
whereas in the previously studied network they are de- 
termined by global quantities. This indeed highlights the 
influence of network topology on the finite- volume correc- 
tions. 

We have carried out a set of numerical tests in or- 
der to verify the accuracy of the theoretical predictions 
for the relative error [as given by Eq. (54) together with 
Eq. (70)] and for the absolute differences between the 
reduced mesoscopic and macroscopic concentrations [as 
given by Eq. (55) together with Eqs. (69) and (70)]. The 
two tests are carried out on a six substrate network and 
are as follows. 

Test 1: Heterogeneous catalytic constants and homo- 
geneous Michaelis-Menten constants. The catalytic con- 
stant is chosen to vary with species i by imposing k!^ = 
0.5 X 1.3^~* while the Michaelis-Menten constants are all 
made equal K\j = 0.1. The other relevant parameter val- 
ues are Mn = 0.002, fci = 10, [E\]V = 2, and V = 200. 
Figure 4(a) shows a comparison between the theory and 
the stochastic simulations. Note that the relative error 
for each species is different, in contrast to that for the 
previous network [inset of Fig. 3(b)]. 

Test 2: Heterogeneous Michaelis-Menten constants 
and homogeneous catalytic constants. The catalytic con- 
stants are all made equal, k\ = 1, while the Michaelis- 
Menten constants are made heterogeneous by choos- 
ing the complex dissociation constants to be substrate- 
specific, k^_i = 9 X 1.3^"*. The other parameters are 
fcin = 0.002, fci = 200, [E\]V = 1 and F = 400. Note 
that in this case the reduced macroscopic concentrations 
arc all equal, [Ai\/K\j = 4, since Eq. (69) only de- 
pends on the catalytic constants. In contrast, the re- 
duced mesoscopic concentrations are predicted and con- 
firmed by simulation [Fig. 4(b)] to vary from one species 
to another. This is due to the explicit dependence on the 
Michaelis-Menten constant in Eq. (70). 



V. DISCUSSION 

In this article we have developed a new approach to the 
calculation of finite-volume corrections to the substrate 
concentrations as predicted by the rate equation theory 
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FIG. 4: Multi-substrate network with multi-enzyme species. 
(a) Heterogeneous catalytic constants and homogeneous 
Michaelis-Menten constants; (b) Heterogeneous Michaelis- 
Menten constants and homogeneous catalytic constants. The- 
oretical predictions are shown as open circles while simulation 
results are the data points. Theory predicts that the relative 
error for this network is generally substrate species-specific 
unlike for the multi-substrate, single-enzyme network where 
it is the same for all species. The insets confirm this via 
simulation. For case (b) theory also predicts that the re- 
duced macroscopic concentrations of all substrate species are 
equal but the reduced mesoscopic concentrations are different 
for each substrate. This is once again confirmed by simula- 
tions which indeed show that the mesoscopic concentrations 
increase as we go further "downstream" in the network. De- 
tailed parameter values can be found in the text. 



of enzyme-catalyzed networks. This theory is comple- 
mentary to the general theory of EMREs recently devel- 
oped by Grima.^ The EMRE theory gives explicit analyt- 
ical results only for relatively simple biochemical circuits; 
we have verified that for one such simple case, i.e., the 
single-substrate, single-enzyme reaction, the present the- 
ory gives the same result as the EMRE. This is an impor- 
tant benchmark for both theories since they were derived 
using completely different methods (the Poisson Repre- 
sentation and the system-size expansion). A noteworthy 
achievement of the present theory over the EMRE the- 
ory is that it can produce analytical results even for very 
complex networks involving arbitrarily large numbers of 
species; this comes about by imposing a separation of 



timescales of all substrate and complex species in the 
network which results in a very substantial dimensional 
reduction of the CME, a result akin to that obtained from 
applying the QSSA on the rate equations. To our knowl- 
edge, this is the first time that the "mesoscopic QSSA" 
has been used to obtain a simple analytical picture of 
the stochastic dynamics of large networks in finite vol- 
umes, i.e., considering effects beyond the conventional 
fiuctuation-dissipation theorem or equivalently beyond 
those which can be predicted using the linear-noise ap- 
proximation. Previous studies employing the mesoscopic 
QSSA have been principally occupied with the need to 
increase the speed of stochastic simulation. ^^"^^ 

In particular, using our approach, we have conclusively 
shown that the predictions of the rate equation models 
for the substrate concentrations will always underesti- 
mate the real substrate concentrations of a network con- 
fined in a finite volume. A generic prediction is that 
the size of the differences increases with decreasing sub- 
cellular volume, decreasing Michaelis-Menten constants 
and increasing enzyme saturation. The topology of the 
network plays a determining role in the magnitude of 
these corrections and how they vary according to the "po- 
sition" of the substrate species in the network. For ex- 
ample for the multi-substrate, single-enzyme network we 
find that the mesoscopic concentrations of all substrate 
species are larger than the RE predictions by the same 
factor (a global effect) whereas in the multi-substrate, 
multi-enzyme network the factor is species specific (a lo- 
cal effect). 

We note that in this study we have not considered 
the effects of diffusion and of finite particle size on 
the stochastic kinetics; for the case of a simple single- 
substrate, single-enzyme reaction, these factors have 
been shown to induce a renormalization of rate constants 
which is apparent even at deterministic length scales. ^^ 
This suggests that the finite- volume corrections that we 
estimated in this paper constitute only a lower-bound on 
the actual differences between the real mesoscopic de- 
scription (which takes into account intrinsic noise and 
noise stemming from both factors) and the conventional 
RE description (which ignores both factors and also in- 
trinsic noise). 

Concluding, we have presented a complete method of 
calculation by which one can study the properties of large 
scale stochastic enzyme kinetic networks. We believe that 
the ability of the method to give simple expressions de- 
scribing the kinetics of such large-scale networks is un- 
precedented and hence may lead to new insight into the 
effects of noise in biologically relevant networks. 
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Appendix A: Derivation of the Poisson 

representation from the moment generating function 

equation 

Consider a birth-death Markov process with moment 
generating function G{z), which evolves as dtG[z) = 
-RPG, where 7^'^ = 'Rp{z,\Iy^. Note that the generic 
form of Tip has ah the derivatives acting to the right. A 
complex Poisson representation can be obtained by ex- 
panding G 

G{z,t)= fdaf{a,t)e^'^'^'^^"^^. (Al) 

Note that here the summation in the exponent extends 
over all components of z and a. The representation is 
carried out in intensive variables a^. Under the condi- 
tion that / is sufficiently compact, / satisfies a Liouvillc 
equation dtf{a) = TZf, where TZ = TZ{a, Vq). The rela- 
tion between TZ and TZ can be clarified by differentiating 
(Al) with respect to t 

dtG{z)^TZ^G{z) = /'dQ:/(a)7^^(2,V;,)e^(^•-l)"'^ 

(A2) 

= f daf{a)TZ^{l + V-^Va,Va)e^'-''-^^"'^. 

(A3) 

If / is sufficiently compact, such that by partial integra- 
tion all boundary terms vanish, we obtain 



dtG{z) ^ I doc (7e^(l - l/-iV„, Fa)/(a))e 



(A4) 



S = —{zCi ^z'^)t,, where the steady-state distribution 
TT is obtained from 



ClTT = 0. 



(Bl) 



The quantity is related to the spectrum of the fast meso- 
scopic variables z described by Ci at zero frequency. It 
enters the reduced evolution as a transport coefficient as 
a consequence of the fast relaxation of the variable z. 
To see this we follow the derivation in^'^ and consider 
{z£^^z'^)jr ~ J dzzjC^^{l — P)z^'K, since Pztt = 0, as 
required. Now consider: 

dsexp{-Cis) = -— exp(-£is)|[^ = C{^{\ - F), 

(B2) 

where we have used that P = lim^-i-oo exp(— £is). Notic- 
ing further that exp(— £is)z-^7r is a solution to — 9^/ = 
L\J with initial condition z^ -k^ we obtain: 



{zL\^ 



z )-K = 



ds{z{s)z'^{Q)),. 



(B3) 



Thus if the adiabatic elimination is carried out in the 
Poisson representation, the transport coefficients are con- 
nected only to the spectrum of the mesoscopic Liouvil- 
lian £i, which is independent of the macroscopic shot 
noise contribution. Adiabatic fluctuations are necessarily 
of Ornstcin-Uhlcnbeck form, as shown in the main text. 
Therefore we consider (Bl) in the form of the potential 
condition: 



{Mz - DV^) TT = 0, 



(B4) 



Therefore the relation between TZ'~^ and TZ is given by 

TZ{a.,ya) = n^{l-V-^\Ia,Va.). (A5) 

Thus formally the Liouville operator TZ can be obtained 
by replacing each Zi by (1 — V~^dl^) and each dl by Vai 
in7^"^. 

Appendix B: Explicit form of the correlators 

The reduced description of enzymatic networks in- 
volves the matrix of transport coefficient, defined by 



from which we can easily verify that {zz^)t^ = 
— (DM^^)"^; the latter is supposed to yield a symmet- 
ric expression. Then using the regression theorem, we 
have: 



POO 

{zC^^z^)^ = ds exp(-Ms)(2z^), 

Jo 



-M^^D^M"'^. 



(B5) 



Hence it follows that S = M ^D^M""^. In a sim- 
ilar manner, one can also prove that {zC^ ^D-n- "= 



^D) 



-M" 
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